The same physics underlying SGRs, AXPs, and radio pulsars 
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Unexpected sign, significant magnitude and time-variation of frequency second derivative exist 
not only in singular radio pulsars but also in Soft Gamma repeaters (SGRs) and Anomalous X-ray 
pulsars (AXPs). This paper shows that the these phenomena are related, and can be interpreted by 
a simple unified model, long-term orbital effect. Thus many of previous "isolated" pulsars may be 
binary pulsars, i.e., orbital period P b w (47, 72)min for AXP IE 2259+586, and P b « (20, 34)min 
for PSR J1614-5047, P b » (3.6, 6.4)min for SGR 1900+14, and P b » (1.5, 5.8)min for SGR 1806-20. 
In addition to X-ray pulsar, IE 1207.4-5209, these two SGRs may be the new ultra-compact binary 
pulsars with orbital period of a few minutes. There might be more binary pulsars that suitable for 
gravitational wave detection than we had thought. 

PACS numbers: 04.30.Db, 97.60. Jd,98.70.Rz 



I. INTRODUCTION 

Pulsars are powered by rotational kinetic energy and 
lose energy by accelerating particle winds and by emit- 
ting electromagnetic radiation at their rotational fre- 
quency, v. The slowdown is usually described by 

v = -kv 11 , (1) 

where k is a positive constant which determined by the 
moment of inertia and the magnetic dipole moment of the 
pulsar and n is the braking index. By Eq.Q we have, 



v jv = nv jv , 



(2) 



n = 3 for constant spin-dipole angle and dipole mo- 
ment. Distortion of the magnetic field lines in the ra- 
dial direction from that of a pure dipole, pulsar wind, 
and time-variable effective magnetic moment results in 

1 < n< 3 HH. 

However a large numbers of first and second pulse 
derivatives of pulsars 3] shows that majority pulsars dif- 
fer from n = 3 significantly, 

\VobslVobs\ > \vjv\ = 3|*>/H ■ 

Most frequency second derivative reported in the litera- 
ture deviate significantly from the simple dipole braking, 
and is usually treated as long-term timing noise. The 
main characteristics of the frequency second derivative 
are: (1) the magnitude of it depends on the length of 
the data span; (2) the sign of it can be both positive 
and negative; (3) the magnitude of i) b s can be orders of 
magnitude larger than that expected by magnetic dipole 
radiation. 

Gong citeGong05a introduced long-term orbital effect 
to explain the puzzles of radio quiet neutron star IE 
1207.4-5209, in which the discrepancy between the mea- 
sured i> bs and the magnetic dipole radiation induced v is 



attributed to the long-term orbital effect. This model ac- 
tually provides a mechanism that can interpret the three 
characteristics of timing noise. 

This paper extend the model in three aspects firstly to 
different kinds of pulsars, SGRs, AXPs and radio pulsars; 
secondly to higher order of frequency derivatives, which 
can explain puzzles in braking index and timing noise, 
and finally the orbital period of 5 different pulsars are 
estimated which can put the model under extensive test. 



II. THE MODEL 

The time for the pulsed light to travel across the pro- 
jection of the orbit into observer's line of sight from the 
instantaneous position of the pulsar is 



r sin i 



sin(w + /) . 



(3) 



where c is speed of light, r is the distance between the fo- 
cus and the pulsar, / is the true anomaly, w is the angular 
distance of periastron from the node, i is orbital inclina- 
tion, as shown in Fig 1. The second of the phenomena 
due to orbital motion is the change of pulse frequency, 



Av 

v 



K[cos(u> + /) + ecosw] 



(4) 



where K = 2ira p sini/[cPf,(l — e 2 ) 1 / 2 ] is the semi- 
amplitude, e, Pb, a p are eccentricity, orbital period, and 
pulsar semi-major axis respectively. 

What if a pulsar is in a binary system, however the ef- 
fect of Roemer time delay and Doppler shift, as given in 
Eq.(J3J) and Eq.Q respectively, has not been measured? 
In such circumstance what the observer measured is nei- 
ther as a true isolated pulsar, no as an usual binary pul- 
sar (which has measured Roemer delay or Doppler shift 
directly) . 

When a pulsar has a companion, the time received by 
the observer (Baryon centric time) is, 
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(5) 
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where t p is the proper time of the pulsar, and z/c is 
dependent of Kepler equation, 



E — e sin E = M = nt . 



(6) 



where M, E and n are mean anomaly, eccentric anomaly 
and mean angular velocity respectively. Notice that t is 
the time of periastron passage, which is uniform. 

Obviously for a true isolated pulsar, we have z/c = 
in Eq.(|5j), thus t\, = t p , which means both tb and t p are 
uniform. But for a binary pulsar system, tb is no longer 
uniform, whereas t p is still uniform. 

Therefore, the proper time of the pulsar, t p , can be 
used to replace the uniform time, t of Eq.©, then we 
have M = nt p . 

If Av of Eq.(0} is averaged over one orbit period by 
the measured time, tb, then it gives 



< Av >= 



1 

Pb 



Avdtb 



1 

Pb 



Av{dt p H — dt p ) 



1 / * " z x 
. Ai'—dl, 



b Jo 



xKv 



PbJo 



Ai/cos(w + E)Edt p 



Pb Jo 



2ji 



[cos(w + /) + ecosw] cos(w + E)dE 



The derivative of Av can be obtained from Eq. (J7J) and 
Eq.@ B, 



*i = /?-(i-0. 

a 



(9) 



where a is the semi major axis of the orbit, and £ = 
2(i+ e 2)(i- e 2/ 4) + t^2- By the expression of a/ a |5|, it is 
related to the advance of precession of periastron, loqr 
by 



a/a _ M X M 2 ! + ■ 



3M 2 (1-e 2 ) 3 / 2 



(2 



3M 2 
2Mi 



y) ' 

(. 10 ) 

where M%, and Mi are the mass of the pulsar, companion, 
respectively and M is the total mass of a binary system. 
P x , P y , Q x , Qy are sine and cosine functions of oj (the 
angle distance of periastron from the node) and (the 
longitude of the ascending node) @. Notice that d/a is 
a long periodic term, whereas, lugr is a secular term. 
Considering Eq.©, the observational v b s is given 



v bs = i> + i>L , 



(11) 



where v is the intrinsic one, which caused by magnetic 
dipole radiation and vl is caused by orbital effect. 



III. APPLICATION TO SGRS 



xKv 



7r(l 



-) + 0(e 4 ). 



(7) 



where x is the projected semi-major axis, x = a p smi/c. 

In practical observation, an observer may average Av 
from to T (T > Pj) through tb, the time received 
by observer, without knowing the orbital period, Pb, at 
all. However if the pulsar measured is truly in a binary 
system, Pb, will affect the averaged result, as given in 
Eq.0, thus the averaged Av given by observer is, 



< Av >- 



T 



Avdtb = -[\ Avdtb + 



P b N 



Pb(N-l) 



Avdh + I Avdh] =0 + o{(3^) , (8) 

NP h 1 



where /? = 7r(l - e 2 / A)xKv / P b . Eq.JTJ and Eq.© indi- 
cate that if a pulsar is in a binary system, then Av (for 
convenience brackets, <, > are not used hereafter) mea- 
sured by the observer is actually contaminated by orbital 
effect, (3. Comparatively, for a truly isolated pulsar, there 
is no orbital effect, and thus f3 = 0. And for a pulsar that 
has already been recognized as in a binary system, the 
effect of p has been absorbed by binary parameters, such 
as, Pb, e, and logr, the well known advance of periastron 
given by General Relativity. 



SGRs and AXPs are believed to be magnetars. Like 
radio pulsars these sources are spinning down and show 
timing noise, however at magnitude on average 100 times 
larger. 

The derivatives of SCR 1900+14 satisfies following re- 
lation, 



j V bs 
Vobs 



(3) 
| V obs | 

Vobs 



(4) 
obs 
(3) I 
obs 



(12) 



The derivatives of SGR 1806-20 satisfies following rela- 
tion, 



(3) 

| — | » |^| 

Vobs ^obs 



obs I 
(3) 1 
y obs 



(13) 



The magnetic dipole radiation satisfies 

v/v = 2,v/v, v^/v^hv/v, v {4) /v {3) = Iv/v. (14) 

Thus it seems that the relation, Ea. (|12|) . can be explained 
by magnetic dipole radiation. However the magnitude of 
v/v is much smaller than observational ratio given by 
Ea.lfT2)l. Therefore, the relation Eq. lQ cannot be ex- 
plained by magnetic dipole radiation, neither is Ea. (|13fl . 
Through Ea. (|ll|l following relation can be obtained, 



Vobs _ 3v v 
Vobs 



: 1~ Wi 

V Vobs Vobs 



(15) 
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where 6j\ = vl/vl- The first term at the right hand side 
of Ea. l(T5|l corresponds to the magnetic dipole radiation 
and the second one corresponds to the orbital effect which 
can change sign and have much larger magnitude than 
that of the first term at right hand side of Eq. (|15l) . This 
actually explains why \n\ 3> 3 is inevitable when the 
second term is ignored. 

Similarly by Ea. Hll|) we can have, 



(3) K • •■ 

V obs _§V V 

Vobs 



V V obs 



■ uj 2 - 



Vobs 



(16) 



where Co 2 



- „( 3 ) 



/i>L- Again the left hand side of Ea.l |16(l 



can be explained by the second term at the right hand 
side which is dominant. lu\ and 6j 2 are generally in order 
of magnitude, Co or Q. However since they are all long pe- 
riodic terms which dependent of ui and f2, at certain time 
the trigonometric function of u) and O may cause a rela- 
tive large or small values in the derivative or derivatives 
of pulse frequency, and results the discrepancy between 
En.ipjl and Eq.JTSJ). 

The magnetic dipole radiation induced frequency first 
derivative is always negative; whereas, the long-term or- 
bital effect can cause both negative and positive vl . The 
fact that most pulsars have negative frequency deriva- 
tive indicates that \v\ > \i>i\. In other words, for most 
pulsars £> bs is dominated by v, thus v bs and v have the 
same sign. Therefore we can assume: v = av Q bs, and in 
turn i*L = (1 — a)P bsi where a > 0.5. For convenience 
define a = 1 — a, notice that a can both be positive and 
negative. Ea. ((T5|) can be rewritten, 



Vobs 
Vobs 



3(1 - afv 



(17) 



Actually a can be obtained from observation, from which 
u)\ of Eci. (|15|) can be obtained. However Coi may deviate 
from Co G r, 



ujgr = 



(18) 



By assuming 7 (typically ~ 10) one can obtain Cogr, from 
which the orbital period can be estimated, 



Pb = 2tt[ 



3(GM) 2 / 3 
c 2 (l - e 2 )uj GR s 



13/5 



(19) 



where G is the gravitational constant. Having p,, the 
semi-major axis of orbit, a of Ea. (|20|l . can be obtained. 
Finally putting Coi, a, p,, as well as estimated M 2 
(Mi = IAMq) into following equation (which is given 
by rewritten Eq.©), 



GMv a 



avobs =vl = 



2i:c 2 a a 



GMv 
2~KC 2 a 



poji 



(20) 



where p = 7rsin i(M 2 /M) 2 (l - e 2 /4)/Vl - e 2 , we can 
obtain p. Then one can adjust the companion mass, M 2 , 
and check whether I sinil < 1 is satisfied or not in the 



expression of p. Moreover putting M 2 , sini and Pf, into 
Eq.© one can obtain Av/v induced by the orbital effect, 
and further check whether it can be consistent with Au/u 
given by observation or not. 

In other words, one can adjust two parameters, 7 and 
M 2 to satisfy the three constraints, a, sini and Av/v. 

By the observation of SGR 1900+14, P obs = 8.2(6) x 
10~ n s s" 1 in May 31-Jun 9, 1998; and P obs = 5.93(3) x 
10" n s s" 1 in Aug 28-Oct 8, 1999 @, v obs varies sig- 
nificantly (v obs = l/P obs ). By a = 0.32, and through 
Ea. (fr7|) and Eq. lfr^|) . we can obtain P b k, 6.4min. By 
assuming Mi = 1.4M©, M 2 = 1.4M 7 = 10, we 
have sini = 0.2 and x = 0.04 through the definition 
of p. Another solution is P& w 3.6min, Mi — 1.4M Q , 
M 2 = 0.5M Q and 7 = 20 as shown in Table 1. 

Nevertheless \a\ < 0.5 guarantees that the magnitude 
of the frequency first derivative caused by the long-term 
orbital effect is smaller than that of the intrinsic spin 
down. This means the long-term orbital effect model 
is not contradictory to the assumption that SGRs have 
high magnetic field^while it is not contradictory to other 
possibilities either [7||8j- 

On the other hand, the second order frequency deriva- 
tives caused by the long-term orbital effect is much larger 
than that of the magnetic dipole radiation induced one 
for SGRs, in other words, v ob s is dominated by Ul. And 
since i>l can change sign at time scale ~ 2tt/Cjgr\ it is ex- 
pected that v obs (also v^ b \, v^ s ) will change sign at time 
scale ~ 2it/ujgr- This can be tested by observation. 



IV. APPLICATION TO AXPS AND RADIO 
PULSARS 

The variation of v b s and i) b s between 1990 Jan- 
uary and 1998 December of 11 pulsars is measured us- 
ing ATNF Parkes radio telescope 9J. The signs of i) b s 
of PSR J1614-5047 and PSR 1341-6220 change for two 
times, which have been attributed to glitch. 

The change on v ohs of PSR J1614-5047 is about 1%, 
thus a — 0.01, and similarly through Ea.l|T5 )l -Ea.l|TS ]l . 
the orbital period of PSR J1614-5047 can be estimated, 
Pb ph 34min. This result corresponds to M 2 — 0.7M Q , 
sini = 0.2 and x = 0.05s as shown in Table 1. 

The measured time scale of change sign on i) bs of PSR 
J1614-5047 is ~ 3.2yr, which is consistent 2tt/logr ~ 3 
yr corresponding to p, « 34min. 

By Eq.© and Eq.©, Av and vl contains orbital ele- 
ments, a, e, and i, which are all long-periodic under the 
Spin-Orbit coupling model. Thus Av and vl are trigono- 
metric functions of ui and (Gong 2005a). In such case ui 
and Q are non-uniform which may induce abrupt change 
in Av and z/£, and therefore mimic glitches in pulsars. 
The absence of glitches in binary pulsars may be due 
to that the additional times delay caused by long-term 
orbital effect can be absorbed by the uncertainties in pa- 
rameters such as P b , Cjgr and e. Whereas, for isolated 
pulsars, the only possible absorbtion of long-term orbital 
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effect is by rotational parameters, v Q b s and v bs, etc. This 
explains why glitches always happens in young isolated 
pulsars. As shown in Table 1, these pulsars usually have 
short orbital period, i.e., from a few minute to 10 min, 
which can make significant timing noise, but the small 
semi-major axis, x ~ 10~ 2 to 1CT 3 , as shown in Table 1, 
prevents them from being observed as binary pulsars di- 
rectly. 

The timing noise parameter is defined as A(t) = 
log(\i>\t 3 )/6v, and A 8 means A(t = 10 s ). By Eq.©, 
the second derivatives of Vl is given, 



GMv 
2irc 2 a 



(Ml ) 



(21) 



Eq. (|21|l indicates that Vl oc p, and by the definition of p, 

we have vl oc sin 2 i, Mf , This explains why PSR 

J1314-6220 has the minimum timing noise, Ag w 0.0012, 
and SGR 1806-20 has the maximum, Ag = 5.5 among 
the five pulsars shown in Table 1. 



DISCUSSION 



2. Magnitude, sign and variation of frequency second 
derivative, u bs- 

3. The relationship of ratios of derivatives of pulse 
frequency like i> obs /v obs , v K ob ' s /v obs . 

4. Why unexpected sign of v b 8 appears much more 
often than that of v a bs- 

5. Why timing noise parameters, Ag, of SGRs are 
much larger than that of radio pulsars. 

There are two predictions from to the new model. v bs, 



,(3) 



and 



obs 



of SGRs, AXPs and radio pulsars should 



change sign at time scale ~ 2t{/Cjqr. The test of this 
prediction may tell us whether timing noise is caused by 
the long term orbital effect or not. 

Another prediction is orbital period of different pulsars 
listed in Table 1. The two SGRs may be ultra-compact 
binary with orbital period of a few minutes, this is a 
natural extension from the binaries with orbital period 
of lO^^min. If confirmed the population of sources for 
gravitational wave detectors, like LIGO and LISA, may 
increase considerably. 



Different timing behaviors shown in different kinds of 
pulsars may caused by the same physics. The long-term 
orbital effect provides possible explanations to following 
phenomena: 

1. Magnitude and time scale on the variation of pulse 
frequency, Av. 
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TABLE I: Calculated rotational parameters and estimated 
orbital parameters of possible binary pulsars 



pulsars 


n 


A 8 




refs 


P 6 (min) m 2 (M ) 


sin i 


X 


7 


e 


IE 2259+586 


3302 


0.42 


0.10 


1 


72 


0.1 


0.6 


0.06 


10 


0.2 












47 


0.1 


0.4 


0.03 


20 


0.2 


SGR 1806-20 


-826 


5.5 


0.32 


2 


5.8 


0.5 


0.4 


0.03 


1 


0.2 a 












1.5 


0.1 


0.4 


0.03 


1 


0.2 6 


SGR 1900+14 


-3921 


4.5 


0.30 


3 


6.4 


1.4 


0.2 


0.04 


10 


0.1 












3.6 


0.5 


0.2 


0.01 


20 


0.1 


PSR J1341-6220 


-1 


0.032 


0.0012 


4 


27 


0.01 


0.2 


0.001 


10 


0.1 












46 


0.02 


0.3 


0.004 


4 


0.1 


PSR J1614-5047 


-81 


0.12 


0.010 


4 


20 


0.1 


0.2 


0.01 


10 


0.1 












34 


0.7 


0.2 


0.05 


5 


0.1 



a and b correspond to observational data of 1999 and 2000 

respectively. Notice that the orbital parameters of SGR 
1806—20 are estimated through u)2 of Eg. 11611 . instead of u>\ 
of Eg. H5H as the rest of pulsars in this table. u>2 is very 
close to lugr, which corresponds to 7 ~ 1. The test of the 
binary nature of these pulsars can be performed by setting 
Voba = v, and v bs = v, in other words, v oba and v b s satisfies 
the expectation of magnetic dipole radiation, and then use 

orbital parameters to fit the quasi-sinusoidal residual. 
Orbital period, Pb, predicted in this table is obtained in the 
case that the mass of the pulsar is Mi = 1.4Mq. The error 

of the estimated Pb mainly comes from the assumption 
ugr = 7<^i, typically the error in 7 is less than 10, thus by 

Eq. ll9lL the error of Pb is P b + f^ 2 (where oT = 3P;, and 
<t2 = 3Pf,/4). rijAs and Vhjvobs are obtained by references: 
1 Kaspi et al 0; 2 Woods et allp; 3 Woods et alEl; 4 
Wang et aldf 



